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ABSTRACT 

We present the direct imaging discovery of an extrasolar planet, or possible 
low-mass brown dwarf, at a projected separation of 55 ± 2 AU (1'.'058 ± 0'/007) 
from the B9-type star k And. The planet was detected with Subaru/HiCIAO 
during the SEEDS survey, and confirmed as a bound companion via common 
proper motion measurements. Observed near-infrared magnitudes of J = 16.3 
± 0.3, H = 15.2 ± 0.2, K s = 14.6 ± 0.4, and L' = 13.12 ± 0.09 indicate a 
temperature of ~1700 K. The galactic kinematics of the host star are consistent 
with membership in the Columba association, implying a corresponding age of 
301^q Myr. The system age, combined with the companion photometry, points 
to a model-dependent companion mass ~12.8 M Jup . The host star's estimated 
mass of 2.4-2.5 M & places it among the most massive stars ever known to harbor 
an extrasolar planet or low-mass brown dwarf. While the mass of the compan- 
ion is close to the deuterium burning limit, its mass ratio, orbital separation, 
and likely planet-like formation scenario imply that it may be best defined as a 
'Super- Jupiter' with properties similar to other recently discovered companions 
to massive stars. 

Subject headings: planets and satellites: detection — stars: massive — brown dwarfs 
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Introduction 



Stellar mass is emerging as one of the most important parameters in determining the 
properties of planetary systems, along with stellar metallicity. Radial velocity surveys have 
i ndicated that the fr equency of giant planets increases with the mass of the stellar host 



( 1 Johnson et al. 



20101 ). and many of the roughly d ozen exoplanets that have been directly 



imaged so far have had A- type stellar hosts (e.g.. 



Marois et al. 



2008 



Lagrange et al. 



20091 ) 



despite such large stars being in the small minority of surveyed targets. Th ese results hav e 



Janson et al. 



motiva ted targeted imaging surveys for planets around massive stars (e.g 
201 lbl ). The increase in planet frequency with host star mass can be readily explained 
theoretically, t hrough the considerat ion that more massive stars are likely to have more 



massive disks ( jMordasini et al. 



20121 ). On the other hand, massive stars also feature an 
increased intensity of high-energy radiation, which may significantly shorten the disk 
lifetime due to photoevaporation, and thus decrease the time window in which giant planets 
are allowed to form. This raises the question whether there is a maximum stellar mass 
above which giant planets are unable to form. 

In this Letter, we report the discovery of a ~12.8 Mj up companion to the ~2.5 M Q 
star k And, the most massive star to host a directly detected companion below or near the 
planetary mass limit. In the following, we describe the acquisition, reduction and analysis 
of the data used for detection, confirmation, and basic characterization of the companion, 
k And b. 



2. Observations and Data Reduction 

Observations of the k And system extended over a period of seven months (January - 
July 2012) and were carried out on Subaru Telescope. JHK images were collected with 
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A0188 flHavano et al 



2010h coupled with HiCIAO (IHodapp et al 



20081 ). L' measurements 



were carried out with AO 188 coupled with the Infrared Camera and Spectrograph (IRCS; 



Tokunaga et al. 



19981 ). Figure [T] displays the multi- wavelength images of the newly 
discovered companion. Table [1] provides a summary of the experimental measurements, as 
well as relevant values from the literature. Figure |2] shows observed astrometric positions 
of k And b as compared with expected motion of an unrelated background star. The 
sub-sections below describe the observations in greater detail. 



2.1. Subaru HiCIAO/AO!88 JHK Imaging 



We first detected k And b using AO 188 coupl ed with HiCIA O on Subaru Telescope 



20091 ) . The observations used 



on January 1, 2012, as part of the SEEDS survey (ITamural 
a 20" x 20" field of view, 9.5 mas pixels, and an opaque 0"6-diameter coronagraphic mask, 
which helped keep the saturation radius < Qf!5. The images were taken in the near infrared 



(if -band, 



.6 Mm), where youn g substellar objects are expected to be bright with thermal 



radiation (IBaraffe et al. 



imaging (AD I; 



Marois et al. 



2003). Pupil tracking was used to enable angular differential 



20061 ) . 



Data reduction of the 46 exposures of 5 s revealed, at 23 a confidence, a pointlike source 
at 1"07 separation. Follow-up observations in J (1.3 /xm; 177 exposures of 10s), H (1.6 /xm; 
171 exposures of 8 s), and K s (2.2 /xm; 135 exposures of 10 s), collected on July 8-9, 2012, 
using the same observing setup, re-detected the source at 6, 28, and 49 a confidence levels, 
respectively. Unsaturated images of the primary, taken immediately before and after each 
filter's observing sequence, and using a neutral density filter (0.866% for H, 1.113% for K s , 
and 0.590% for J) provided photometric calibration. 



To optimize the ADI technique, we first reduced the data using a locally optimized 
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combination of images algorithm (LOCI; lLafreniere et al.l 120071 ). HiCIAO observations of 
M5, combined with distortion-corrected images obtained with the Advanced Camera for 
Surveys (ACS) on the Hubble Space Telescope, enabled accurat e pixel scale calibrat ion 



to within 0.2%; the ACS astrometric calibration was based on 



van der Marel et al. 



( 120071 ). Figured] (left and middle) presents a JHK false-color image and corresponding 



signal-to- noise (S/N) map after the ADI/LOCI data reduction. 

Given the relatively high S/N ratios and the known difficulties in quantifying the 



impact of 
reduction 



OCI on planet ph otometry and astrometry, we also performed a classical ADI 



Marois et al 



20061 ) with mean-based point-spread function (PSF) estimation 
and frame co-adding. Unsharp masking on the spatial scale of 35 pixels (~ 7FWHM) 
was applied to the final image to flatten the residual background. The planet signal was 
recovered with S/N ratios comparable (within 10%) to the LOCI reduction for all the July 
data sets. For the somewhat lower quality January data, the measured S/N reduced from 
about 23 a to 7 a. 

To achieve unbiased photometry and astrometry, we extracted the combined k And 
PSF (S/N > 1000) from the neutral density images, and placed it on an empty image 
frame at the location of k And b. Applying the same unsharp masking and ADI reduction 
to this data as we did for the science data, we simulated the parallactic angle evolution, 
as recorded in the science frames. The resulting processed PSF acted as the photometric 
and astrometric reference for k And b. The only non-linear step in this process was the 
median-based unsharp masking, but the large spatial scale (~ 7FWHM) ensured that 
subtraction effects were minimal. 

We calibrated the astrometry by cross-correlating the k And b signal with the 



processed calibration PSF . We estimatec 



FWHM/{S/N), following 



Cameron et al 



;he un certainty in the k And b center to be 



( 120081 ). The uncertainties in the final relative 
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astrometry were dominated by our ability to determine the host star center, which was 
carried out through Moffat fitting of each individual exposure. We conservatively estimated 
the uncertainty of the Moffat fit at 0.75 pixels (7 mas). For confirmation, we applied Moffat 
fitting and peak fitting to unsaturated data of k And and found that the methods agreed 
at the 0.5 a level. The photometric uncertainties were calculated as a combination of (1) 
representative noise in an annulus, centered on the host star, with a radius equal to the 
companion, (2) photometric variability in the neutral density calibration images, which 
yielded effective accuracies of 7-11% for the combined datasets, and (3) uncertainties in the 
JHK magnitudes of k And. 



2.2. Subaru IRCS/A0188 V Imaging 



On July 28, 2012, we followed the JHK observations with L'-band observations (3.8 /xm; 
50 exposures of 30 s) using A0188 coupled with the Infrared Camera and Spectrograph on 
Subaru Telescope. We employed a 10'.'5 x lO^ field of view, 20.6 mas pixel scale and no 
coronagraph. The host star saturated out to ~ O'.'l. The dithered observations, carried out 
in ADI mode, were divided into two identical sequence s bracketing observ ations of the star 



HR 8799, which provided the photometric calibration ( iMarois et al. 



20081 ). Observations 



of a third star, S 810-A, were collecte d before the science observations as a secondary 



calibration check (ILeggett et al. 



2003|). 



We sky-subtracted each image using a median combination of frames taken at the 
other dither positions. To help maximize the hig h-contrast sensitivi ty, we processed the 



data using an "adaptive" LOCI process (A-LOCI; 



Currie et al 



20121 ). We also employed a 



moving pixel mask, where the LOCI algorith m is prevented from usi ng, in PSF construction 



pixels lying within the subtraction zone (see 
(right) shows the final image. 



Lafreniere et al. 



20071 for details). Figure Q] 
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To quantify the k And b throughput, we used fake point sources added to the image and 
processed with the same algorithm settings. As an additional check on our flux calibration, 
we determined the relative brightness between the HR 8799 bed planets (all detected 
at S/N > 7-10) using ide ntical procedures, and confirmed its agreement with published 



values ( jCurrie et al 



20111 ) . The independent calibrations all yielded self-consistent results, 
ensuring confidence in the 22 a detection of k And b in L' . As a final check, we re-processed 
the L'-band data using a more classical ADI method, similar to that described for the JHK 
data set, and achieved consistent results. While the July V astrometry was consistent with 
the July JHK results, we refrained from including it in our proper motion analysis, due to 



our possession of poorer-quality astrometric calibration. 



3. Host Star Properties 



k And is a B9 IV star flWu et al 



201 if ) located at a distance of 52.0 pc ( IPerryman et al. 



19971 ). iFitzpatrick fc Massal ( 120051 ) report a temperature of 11,400 ± 100 K w ith a sub-s olar 



met a 



licity of [Fe/H] = —0.36 ± 0.09, while independent measurements by 



Wu et al. 



( 1201 lh report values of 10,700 ± 300 K and -0.32 ± 0.15. Given the star's spectral 



classification, the measured low metallicity is likely due to the details of the star's accretio n 



and a tmospheric physics, as opposed to a true, initial, low metallicity (IGray Corbally 



20021 ). We e stimate a mass of 2.4- 2.5 M & using the published temperature and evolutionary 



tracks from 



Ekstrom et al. 



( 120121 ). Table [U summarizes the host star properties. 



IJ L 

Zuckerman et al.l ( 120111 ) proposed k And to be a member of the ~30 Myr old Columba 



association. To further investigate k And's likely membership in Columba we: (1) 
i ndependently ca 



( IPerryman et al 



culated its Galactic kinematics from astrometry available in the literature 



199 



Zuckerman et al. 



associations reported in 



Torres et al 



20111 ) and compared these to the young local 



(120081 ) . and (2) calculated its membership probability 
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in these associations using the Bayesian methods of iMalo et al.l ( 120121 ) . Our analyses 
showed that the star's kinematics imply a >95% probability of the star being part of the 
Columba association. 



As an additional check, we compared the k And-B — V color and absolute V magnitude 



(jPerryman et al. 



1997 



van Leeuwen 



20091 ) with members of clusters and associations with 



ages ranging from ~15-700 Myr. These include Lower Ce ntaurus Crux, a Per, Pleiades, 



Com a Ber, Hyades, Praesepe, and young local associations ([Torres et al 



2008 



van Leeuwen 



20091 1 . The color-magnitude analysis showed that k And is consistent with other early-type 



stars having ages ~20 -120 Myr. The results o f our analyses are consistent with the 



conclusions reported in 



Zuckerman et al. 



(120111 ); k And's age range and kinematics suggest 
it is a me mber in the Columb a association. We therefore adopt a system age of 30^0 Myr 



(following 



Marois et al 



20101 1 for all subsequent analyses. 



4. Results 
4.1. Proper Motion Analysis 

Located 52.0 pc from the Sun, k And exhibits proper motion of 83.5 mas/yr 



(jPerryman et al 



19971 ). enabling an effective test to distinguish bound companions from 
unrelated background stars. The k And proper and parallactic motion translate to 76 
mas (~8 HiCIAO pixels) of net movement over the 6 month period between epochs. As 
shown in Figure [2J the companion exhibits common proper motion with the host star, and 
deviates from expected background star motion by 7 a. In addition to this 7 a deviation in 
the magnitude of motion, the observed direction of motion and scatter in astrometry are 
completely inconsistent with that of a background star. 
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4.2. Physical Properties of k And b 



Figure [3] shows that the k And b colors are most consistent with cloudy L dwarfs 
and overlap with several other benchmark exoplanets and low-mass companions, including 
HR 8799 bed, AB Pic b, and 1RXS1609 b. Figure [j] co mpares k And b co l ors and absolut e 



magn itudes with DUSTY and COND evolutionary tracks (IBaraffe et al. 



2003 



Chabrier et al.l 



20001 ) . as well as low-mass companions around HR 8799 and AB Pic. The plots show 
k And b as well situated between HR 8799 cde and AB Pic b. Its infrared colors are 
slightly bluer tha n those of typical field L dwarfs, possibly indicating a low surface gravity 



( ICruz et al 



20091 ). However, improved photometry is required to confirm whether this color 



deviation is real. 



The estimated temperature of k And b suggests that its atmosph eric properties shoul d 



Chabrier et al. 



200 



align more closely with those of the DUSTY models (see discussions in 
In deriving a mass estimate from this track, we rely on the July if -band magnitude because 

(1) alternative J-band and January ii-band measurements have higher uncertainties, 

(2) ifs-band mass estimates are more sensitive to atmospheric composition (see e.g. 



Janson et al. 



201ll ). and (3) L'-band mass estimates have been less thoroughly tested with 



experimen tal data, and are mo re sensitive to age uncertainties for this age and magnitude 



range (see 



Chabrier et al. 



20001). 



Based on the July if -band magnitude of 



a parallax of 19.2 ± 0.7 (jPerryman et al. 



5.2 ± 0.2, the estimated age of 30^ M J r , 



19971 ). and the DUSTY evolutionary models, we 
calculate a mass of 12.8^q Mj up and a temperature of 1680^20 K. As a consistency check, 
we calculate the predicted JK S L' magnitudes based on the estimated 1680 K temperature, 
the 20-50 Myr system age, and the DUSTY evolutionary models. This yields J = 16.5-16.8, 
K s = 14.2-14.4, and V = 13.1-13.2 mags, all of which are in agreement with our measured 
multiband photometry. Additionally, the two epochs of ii-band photometry are in 
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agreement with one another. Table [T] summarizes the complete properties of k And b. 

While the DUSTY models are likely the more relevant, we estimate a possible 
alternative mass using the COND evolutionary tracks. In this scenario, we determine a 
mass of M ^ !P and a temperature of I640I20 K. More recent evolutionary models by 



Spiegel fc Burrows! ( 120121 ) offer alternative "Warm Start" scenarios that consider formation 
with lower levels of initial entropy. While these models do not consider combinations of 
mass and temperature similar to that of n And b, they do predict generally higher masses 
than that of the DUSTY and COND models. In the case of k And b, such models place 
the most probable mass at a value above the typical deuterium burning limit. While we 
currently adopt a nominal mass estimate of 12. 8^ Mj up for the analyses in this discovery 
paper (based on the DUSTY models), we defer a deeper investigation of companion mass 
for a follow-up paper, where we will focus on a more thorough comparison of multiband 
photometry with synthetic spectra. 



4.3. Orbital Properties of k And b 



We estimate the semimajor axis of k And b from its observed separatio n. Assuming a 



unifor m eccentricity distribution of < e < 1, and random viewing angles, 



Dupuy et al 



( I2010l ) compute a median correction factor between projected separation and semimajor 



axis of Lllg gg. Using this relation, we derive a semimajor axis of 6ll2o AU based on its 
projected separation of 55.2 AU (1'.'07) in January 2012. 



4.4. Possible Secondary Companions 



The if -band sensitivity levels (see Section 2.1) allow us to rule out secondary 
companions with temperatures similar or warmer than that of k And b, for separations 



13 



greater than 0'.'9 (46 projected AU). For the k And b separation (l'.'l) and beyond, we may 
rule out secondary companions with masses > 11.7 Mj up , assuming a 30 Myr system age 
and the DUSTY evolutionary models. 



5. Discussion 



k And is the most massive star to host a directly imaged planet, or brown dwarf 
near the deuterium burning boundary. The mass ratio between k And b and its host is 



0.5 %, similar to the ~ 0.4% ratios of the Pic and HR 8799 planets (ILagrange et al. 



2009 



Marois et al. 



20081 ). In comparison, this value is noticeably smaller t han those of 



reported direct 



(jChauvin et al. 



y ima ged planets around 1RXS 1609 ( jLafreniere et al. 



2008|) and 2M 1207 



20041 ). The projected separation of « And b is also intermediate between 



the two outer planets in HR 8799. The similarities between k And b, Pic, and HR 8799 
could imply a similar formation mechanism, which may be distinct from recently discovered 



brown dw arf companions o: 



GJ 758 B; 



: appr oximately an order of magnitude larger mass ratios (e.g. 



Thalmann et al. 



20091 ) or semimajor axes (e.g., HIP 78530 B; 



Lafreniere et al. 



20 111 ) . Stren gthening the 



models (e.g. 



Rafikov 



possibility of a planet-like formation for Kap And b, theoretical 



20 111 ) show that, for a minimum mass solar nebula, the region of the 
primordial disk where core accretion formation of giant planets can occur overlaps with the 
separation range of k And b. Furthermore, this formation mechanism may be significantly 
enhanced for a star as massive as k And, assuming it had a correspondingly more massive 
protoplanetary disk. Further studies will be needed to more stringently constrain the 
population properties of planets and brown dwarfs on intermediate and wide orbits. 

The best-fit mass of k And b lies just below the deuterium burning limit according 
to conventional evolutionary models, but ma y be above this limit if i nitial entropy at 



formation is lower than such models assume (jSpiegel fc Burrows 



20121 ). This leads to an 
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ambiguity in whether the companion can be classified as an "exoplanet" by the present 
IAU definition. Such a classification scheme can however be misleading, given that k And b 
may well have formed in the same way as previously imaged planets, regardless of whether 
its mass falls just below or above this limit. Indeed, radial velocity studies have shown that 
massive stars tend to have massive pla nets, sometimes with companions having masses 



above the deuterium burning limit (e.g. 



Lovis fc Mayor 



2007) and which appar ently form 



Hekker et al. 



2008|). On the 



a high-mass tail of a lower-mass planetary population (e.g. 
other hand, formation history can be difficult to assess in individual cases. In order to 
avoid these uncertainties, we simply classify k And b as a 'Super- Jupiter', which we take to 
mean a group of objects that includes the previously imaged planets around HR 8799 and 
(3 Pic as well as the most massive radial velocity planets, and which one might suspect have 
formed in a similar way to lower-mass exoplanets, but for which this has not necessarily 
been unambiguously demonstrated. This suggested class includes substellar objects with 
masses at or moderately above the deuterium burning limit, but excludes objects with 
orbital separations well beyond a typical disk truncation radius, or systems with mass ratios 
more indicative of a binary-like formation. 
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(companion) 




Fig. 1. — Left: JHK false-color image of k And b after LOCI/ADI data reduction, for the 
2012 July observations. Center: A corresponding signal-to- noise map created from the left 
frame. The S/N ratio is calculated in concentric annuli around the star. The white plus sign 
in each panel marks the location of the host star k And; the black disks designate the regions 
where field rotation is insufficient for ADI. White features indicate where the signal is roughly 
equally strong in all wavelengths; colored features indicate where the signal is mismatched 
between wavelengths, and is often indicative of residual noise. The lobes around k And b 
result from the Airy pattern produced by the Subaru A0188 system. Right: L'-band image 
of k And b from the 2012 July observations. 
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Fig. 2. — Proper-motion analysis of k And b. The dotted curve designates the predicted par- 
allactic and proper motion between epochs, if the detected January source were a background 
star. The dashed line indicates an example bound, orbital path of k And b consistent with 
the observational data. The diamond symbols represent the predicted January and July as- 
trometric measurements for n And b, if it follows the dashed orbital path, k And b is clearly 
inconsistent with background behavior and instead demonstrates common proper motion 
with the host star. 
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Fig. 3. — Position of k And b colors (red points) with resp ect to reference obje cts. Top plot 



i ncludes benchmark substellar c ompanions: HR 8799 bed (IMarois et al. 



( jChauvin et al. 



2004), AB Pic b ( 



Chauvin et al. 



2005h. 1RXS1609 b (ILafreniere et al. 



2008) , 2M1207 A b 



2008|), 



CD-35 2722 B f|Wahhaj et al 



( iBejar et al. 



face gravity (I Cruz et al 



20111), GSC 06214 b fllreland et al. 



20111), and USCO 108 AB 



2008) . It also con t ains L dwarfs with s pectral features indicative of reduced sur- 



2009 



Faherty et al.ll2012[ ). and Pleiades M-L dwar fs ( IBihain et al. 



2010h. Bottom plo t includes M, 



( jCurrie et al. 



2011 



Skemer et al. 



, and T field dwarfs ( 



20121 ). and 2M1207 b (jChauvin et al. 



>eggett et al. 



2002), HR 8799 bede 



20041 ) . 
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Fig. 4. — k And b colors and absolute magnitudes (red points') co mpared with DUSTY and 
COND evolutionary tracks ( IBaraffe et al.lbooa Ichabrier et al.lbooct ). 
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Table 1. Properties of the k And System 



Property 



Primary 



Companion 



Mass 
T e g 

Spectral Type 
Age (Myr) 
Parallax (mas) 
Fe/H 

logs 

J (mag) 
H (mag) 
K s (mag) 
L' (mag) 
AJ (mag) 
A_ff (mag) 
A_ftT s (mag) 
Mj (mag) 
M H (mag) 
M Ks (mag) 
M L , (mag) 

Astrometry on 1 January 2012 (_H"-band): 
-Proj. scp. (") 

— Proj. sep. (AU) 

— Position angle (°) 
Astrometry on 8 July 2012 (ff-band): 
-Proj. scp. (") 

— Proj. sep. (AU) 

— Position angle (°) 



2.4-2.5 M" a 
11,400 ± 100 K c 
10,700 ± 300 K d 
B9 IV d 

o n +20f 

au -io 
19.2 ± 0.7 g 

-0.36 ± 0.09 c 

-0.32 ± 0.15 d 

4.10 ± 0.03 c 

3.87 ± 0.13 d 

4.6 ± 0.3 h 

4.6 ± 0.2 h 

4.6 ± 0.4 h 



1.0 ± 0.3 1 
1.0 ± 0.2 1 
1.0 ± 0.4 1 



12.81*;° M. Jup b 
1,680^20 Rb 



L2-L8 e 



16.3 ± 0.3 
15.2 ± 0.2 
14.6 ± 0.4 
13.12 ± 0.09 

11.6 ± 0.2 
10.64 ± 0.12 
10.0 ± 0.08 

12.7 ± 0.3 
11.7 ± 0.2 
11.0 ± 0.4 

9.54 ± 0.09 

1.070 ± 0.010 
56 ± 2-i 
55.7 ± 0.6 

1.058 ± 0.007 
55 ± 2' 
56.0 ± 0.4 
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Note. — Photometric values represent Subaru July 2012 measurements, unless 
noted otherwise. 



Calculated usin g the published tempe rature from 



tionary tracks from 



Wu et al 



Ekstrom et al 



I 20111) and evolu- 



1 2013) . 



b Calculated u sing the Jf-band magn itude, estimated re And age, and evolution- 



ary models from 



Chabrier et al 



(2000). 



Fitzpatrick & Massal (2005) 



c Wu et al 



l 20111) 



°Based on measured colors and 



Cruz et al 



I 2009 ) spectral identifications. 



*Zuckerman et al 



( 20111) and 



Marois et al 



(2010) 



g Hipparcos 
h 2MASS 



Perrvman et al 



(1997) 



Skrutskie et al 



(2006) 



'Calculated by the authors, using 2MASS photometry and Hipparcos parallax. 
J Uncertainty is dominated by the host star parallax measurement. 



